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The effects of Umklapp scattering on electronic states are studied in one spatial dimension at 
absolute zero. The model is basically the Hubbard model, where parameters characterizing the 
normal {U) and Umklapp {V) scattering are treated independently. The density of states is 
calculated in the t-matrix approximation by taking only the forward and Umklapp scattering 
into account. It is found that the Umklapp scattering causes the global splitting of the density 
of states. In the presence of sufficiently strong Umklapp scattering, a pole in the t-matrix 
appears in the upper half plane, signalling an instability towards the 'G/2— pairing' ordered 
state (G is the reciprocal lattice vector), whose consequences are studied in the mean field 
approximation. It turns out that this ordered state coexists with spin-density-wave state and 
also brings about Cooper-pairs. A phase diagram is determined in the plane of V and electron 
filling n. 

KEYWORDS: modified t-matrix approximation, Umklapp scattering, density of states, mean field approxima- 
tion, G/2-pairing, G/2-SDW. 



§1. Introduction 

The problems of the metal-insulator transition due to the short-range Coulomb repulsion in the 
lattice system (known as the Mott transition^ ) have attracted much interest. The Mott transition 
is often examined based on the Hubbard mo deli with the hopping integral t and the on-site 
Coulomb repulsion U . At least qualitatively, this phenomenon can be understood: If U /t is large, 
electrons are kept apart each other. Especially for a half-filled band (one electron per one lattice 
site), they get completely localized because it requires large energy comparable with U for them to 
move to another lattice site. As a result, the charge excitation spectrum has a gap known as the 
Mott-Hubbard gap and the system becomes an insulator. On the other hand if [//t is small, it is 
energetically favorable for electrons to spread into the energy band and the system behaves as a 
metal. Hence it is understood that both the Coulomb repulsion and the existence of discrete lattice 
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are the keys to the Mott transition. It was Hubbard who first succeeded in describing the quahtative 
feature of such behaviorsB. By using the Hubbard model, he showed that the Mott-Hubbard gap 
begins to appear as U/t gets larger in the half-filled case. 

If we consider the on-site Coulomb interaction U term in the wave vector representation, we 
have not only normal processes (which conserve the total momentum of two scattering electrons) 
but also Umklapp processes. Since the latter processes originate from the periodicity of the lattice 
and stand for the large-momentum-transfer processes, they play an important role in the Mott 
transition. 

In the space dimension d = I, the exact solutions of the Hubbard model have been obtained by 
using Bethe Ansatz methodB. It was proved that the charge excitation spectrum has a gap and 
the system becomes an insulator for a half-filled band for any strength of Coulomb interaction. 
In this method, however, it is not easy to obtain the physically clear-cut understanding of the 
origin of the Mott-Hubbard gap. On the other hand, it was shown by the bosonization schemei^ 
that the Umklapp scattering produces this gap in a half-filled bandHll'. Moreover, characteristic 
features of the charge excitations in one-dimensional doped Mott insulators, especially the crossover 
between the metallic Drude type of excitation and the Mott-Hubbard gap, have been clarified@) by 
mapping the model first to the g-ology modefli'i^, then to the phase Hamiltoniani^ and finally 
to the massive Thirring model. It has also been disclosed that the Mott-Hubbard gap and other 
associated anomalies (the divergence of the compressibility near half filling, etc.) originate from 
Umklapp processes and how the system begins to show anomalies as the system approaches half 
fillinj'&0'§). 

The emergence of the Mott insulator has also been clarified in the case of (i = oo and it was also 
shown that the density of states has a gap together with the existence of quasi-particle peak before 
the onset of Mott transition0\ In d = oo, however, we cannot distinguish Umklapp processes from 
normal ones because of the irrelevance of the momentum conservation. 

From this point of view, it will be important to clarify the roles of the Umklapp scattering 
especially near half filling. In this paper, we examine the effects of Umklapp processes in detail at 
absolute zero for the simplest d = 1 case, bearing the application to the 2D case in mind. 

In §2, we introduce the model Hamiltonian. In §3, we study the density of states based on the 
perturbation theory with a special emphasis on both the forward and Umklapp scattering in the 
t-matrix approximation. Since we found that the t-matrix has a pole in the upper half plane if 
the Umklapp term becomes dominant, we study the nature of this instability in the mean field 
approximation in §4. Another type of instability associated with Umklapp processes is studied 
here. §5 is devoted to conclusion and discussion. 

In the present paper, we take unit of h = ks = 1- 



2 



§2. Model Hamiltonian 

We start with the one-dimensional Hubbard model with nearest-neighbor hopping: 

H = -t J2 {cLcjV + (h.c.)} + UY^rii^nn - /x^cJ^q^, (2.1) 

<ij>u i i 

where t{> 0) is the transfer integral, Cia{cl^) is the annihilation (creation) operator for the electron 
on the i-th site with spin a, U (> 0) is the Coulomb repulsion and /x is the chemical potential. 
We can classify normal and Umklapp processes by the Fourier transformation, 



H = Ho + W\ (2.2a) 

^0= E T^^kcl^ka, (2.2b) 

-G/2<k<G/2 

ik = ^k- 1^, (2.2c) 

ejk = -2icosA;, (2.2d) 

jjint ^ ^int ^ ^int^ (2.2e) 



Tjint 

-Hat 



Q<q<G I q-G/2<k<G/2 
q-G/2<k'<G/2 



E 



k'+q-Gi 



C-k+q-GlCk]{^-?f) 



-G/2<k<q-G/2 
-G/2<k'<q-G/2 



Tjint 



V_ 
N 



E 

0<q<G 



E 



Cfc'T'^-fc'+g-Gi'^-fe+9j.Cfct + (/i.e.) 



q-G/2<k<G/2 
-G/2<k'<q-G/2 



(2.2g) 



where N, Ri, G and g (or g — G) are the number of lattice sites, the lattice site vector, the 
reciprocal lattice vector equal to 27r and the total momentum of two electrons, respectively. Here 
the lattice constant is taken as unity. Each wave vector k of the electron lies in the first Brillouin 
zone, —G/2 < k < G/2. The cosine band is shown in Fig. 1. The bandwidth W is equal to 4t. 

stands for normal processes which conserve the total momentum q of two electrons, and Hjj* 
stands for Umklapp processes in which it differs by G. Since we are mainly interested in Umklapp 
processes, we set the magnitude of interaction in Hjj^ as an arbitrary parameter V{> 0). In this 
way, we examine the effects of normal and Umklapp scattering separately. 
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Fig. 1. Dispersion relation. Full and thin line show the cosine band and the linearized one, respectively. 

§3. The Density of States in the Modified T-matrix Approximation 

In this section, we study the effects of the Umklapp scattering on the density of states based 
on tlie perturbation theory, at absolute zero T = 0. For simpUcity, we approximate the cosine 
dispersion with the hniarized one: 

ek = v{\k\ - G/4), (3.1a) 
^k = v{\k\-kF), (3.1b) 

where kp is the Fermi wave vector and we put v = 8t/G so that the bandwidth W is left unchanged. 
In the half-filled case, kp is equal to G/4. This dispersion is shown in Fig. 1. 

We separate the dispersion into the branch with k > (right-moving) and with A; < (left- 
moving) and treat only the scattering between two electrons on the same branch - the forward and 
Umklapp scattering, as shown in Fig. 2. However, we neglect backward scattering processes in this 
paper, where the total momentum of two incident electrons is small. 




Fig. 2. (a) Forward and (b) Umklapp scattering processes. Full and dotted arrow lines stand for right-moving 
electron and left-moving electron, respectively. 
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3.1 Particle- Particle Correlation Function and the T-Matrix 

First of all we define the particle-particle correlation function K{q, iuoi) with the total momentum 
q and energy iuji shown in Fig. 3, 

K{q,iuii) = T^^Go{k,ien)Go{-k + q,-ien+iuJi), (3.2a) 

= ylM^ItL^^ (3.2b) 

Go{k,ien) = -. — ^— — , (3.2c) 

«en - Kk 

where Go is the thermal Green function of Hq, loi = 21ttT and e„ = (2n + 1)itT are the Matsubara 
frequencies and /(^) is the Fermi distribution function. Since we consider only the interaction 




K(q,icoi) 
-k+q,-iejj+i«)j 



Fig. 3. Particle-particle correlation function, /^(g, ioji). 



between two electrons on the same branch, + is independent of k and we can calculate 

K{q,iuJi) analytically. Especially, for \q\ ~ 2kF at absolute zero T = 0, 

where 

Xg = v{\q\-2kF). (3.4) 

We note that K{q, z) has one pole on the real axis z = Xq and for \q\ ~ Ikp the residue is small. 
Next we define the t-matrix T{q,iu;i) without Umklapp processes (i.e., V = 0), 

T{q,iuJi) = - ^ . (3.5) 
1 + UK{q,iuJi) 



This is the particle-particle ladder diagrams shown in Fig. 4. By use of K{q,z), eq. ( p.3[ ), we can 
calculate the t-matrix analytically. Especially for \q\ ~ 2kF, 

T{q,z) = -U-Ua ..""^ . , (3.6) 

z - [1 + a)Xq 

where 

- U - U 

"=2^ = 2W ('-^^ 
We note that T{q,z) also has one pole on the real axis and for |g| ~ 2kF the residue is small. 
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Fig. 4. The t-matrix, T{q,iLUi). 



3.2 The Modified T-matrix 

We consider the particle-particle ladder diagram including Umklapp processes (i.e., F 7^ 0). In 
this case the total momentum q of incident two electrons can be changed into q — G{q>0) after the 
Umklapp scattering, i.e., the Umklapp interaction introduces the Bragg reflection of two incident 
electrons as a whole. In this case, the t-matrix is given as follows: 



f = -U- Ukf, (3.8a) 

(3.8b) 



Tn{q,z) Ti2{q,z) 
T2i{q,z) T22{q,z) 



U V 

U= \ ^ ^1, (3.8c) 

. K{q,z) , , ^ 

K=\ ' , 3.8d 

' K(g-G,^) ' 

where the index 1 stands for the total momentum g > and 2 for g — G < 0. We call this matrix 

as a modified t-matrix. This is the particle-particle ladder diagram shown in Fig. 5. Especially for 

q ~ 2kF, Tii{q, z) is given by 

"^^^^^'^^ - ^^^) ' ^'-'"^ 

x±{q) = ^ {(1 + + ± ^(1 + aY{xq - Xq^c? + 462xgXg_G| , (3.9b) 

6^ V V 



and has two poles z = x±{q). In the limit F ^ 0, Tii{q,z) reduces to eq. ( |3.6| ). 
In the half-filled case (G = 4/ci?), Xq-c = —Xq and then 



x±{q) = ±axq, a = J{l+af -h^, (3.10) 
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Fig. 5. The modified t-matrix, T(g, ioij). 



and therefore Tii(g, z) for g ~ 2kF is written as 

' {a — l){a + c) ajg 



2a 



ax„ 



{a + l){a — c) Xq 



2a 



(3.11) 



The renormalization factor a of the velocity v becomes pure imaginary if 6 > 1+a, i.e., V > U+2W. 
Therefore a pole appears in the upper half plane in this case. This fact implies that if the Umklapp 
interaction is sufficiently strong, the system becomes unstable. We will study this instability in the 
next section. In this section, a is assumed to be real. 

3.3 The Density of States 

In this section, we calculate the density of states. First of all, the single-particle spectral weight, 
p{k, e), is given by 



p{k,e) 



—lmG^{k,e) 

IT 



(3.12a) 

= pcik,e) + pp{k,e), (3.12b) 

where G is the full Green function with the proper self energy S, G^{k, e) = [e — — S^(A;, e)]~^. 
The superscript R stands for the retarded function. In eq. (3.12), pc{k, e) stands for the contribution 
from continuum part with ImS^(fc, e) ^ 0, 

1 hnE^''{k,c) 



TT {e - a - ReS«(fe, e)}2 + {ImE«(A;, e)}^ ' 
and Pp{k, e) for the contribution from the poles with ImE^(/u, e) = 0, 

Pp{k,e) = 6{e-Ck-ReJ:'^{k,e)), 



(3.13) 
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(3.14) 



where 



k=ki{e) 



(3.15) 



and ki{e) (i= 0, 1, 2, • • •) is the solution of e — — ReS^(A;, e) = satisfying ImS^(/c, e) = for 
given e. The summation is taken over all ki{e). 
By use of p{k, e), we obtain the density of states p{e) per lattice site and spin, 



pie) = pcie)+pp{e), 



where 



1 



(3.16) 



(3.17) 



and V = c,p. In our approximation where the scattering processes between two electrons on the 
same branch are taken into account with the linear dispersion eq. (3.1), we can calculate S^(A:,e) 
analytically by use of the t-matrix in both cases of ^ = and V ^ 0. 
3.3.1 The 2nd Order Perturbation in U 

First of all, we analyze the results in the 2nd order perturbation in U but V = In this case, 
the self energy can be written with K{q, z) as follows: 



S(/c, ien) = -TU'^ Go{-k + q, -ien + iLdi)K{q, iui) 

iu)i q 



(3.18) 



After straightforward calculations, we obtain for A; > 0, 



>^ I < < 2kF, 

2 e-(,k + id 



p2 



S^(fc,e) 



c2 _ e2 
\ Jl^k 'ik-kf 



(3.19) 



^a"— — "-x '^f^F < k < G/2. 
Therefore it is seen that the Green function has two poles at e = E±{k), 



E±{k) 



72^ 



(1 ± ^)^^ 



for < < kp, 
for kp < k < 2kF, 



(3.20) 



± a^kpik - 3/cf/2) for 2kF < k < G/2. 



These energy bands for k > are shown in Fig. 6. In the half- filled case, 2kF = G/2 and it is seen 
from eq. ( |3.2[l| ) that these dispersions are linear in k over the first Brillouin zone. 

The density of states consists only of the contribution from the poles, Pp{k, e), which is shown in 
Fig. 7 for the half-filled case with U/W = 2.0. 
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Fig. 6. The quasiparticle energy bands relative to the renormahzed chemical potential in the 2nd order perturbation 
in f7 for A; > and less-than-half-flUed case. The states with e < are occupied. 
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1 
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Fig. 7. The density of states in the second order perturbation in U for the case of half filling and U/W = 2.0. 



3.3.2 The T-matrix Approximation for U 

Next we study the density of states in the t-matrix approximation for U in the absence of V. In 
this case, the self energy can be written with T{q, z) as follows: 



E(A;, ien) = X! X! G^^-k + g, -ie„ + iu>i)T{q, iui). 

iwi 9 



Especially for < /c < 2kF, T,^(k, e) is given by 

ReS^(A;,e) = -aCfc-(e-efe)ln 



e-(l + a)efe 



e-6 



ImS^(A;, e) = -7r(e - 6)sgn(6) for 1 < e/^k < 1 + a. 



(3.21) 

(3.22a) 
(3.22b) 
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Here the contribution from the first term of the r.h.s. of is neglected, because this can be 
included in the chemical potential shift. It is seen that the Green function G^{k,e) has two poles 
at e = E±{k), 

c^£,k for < < kp, 



E±{k) 



where c± 



c±^k for kp < k < 2kp, 

c+ = 1 + 1.19a, 
c_ = 1 - 0.466a. 



(3.23) 



(3.24a) 
(3.24b) 



are determined numerically. There are two quasiparticle energy bands corresponding to these poles, 
which are shown for < k < 2kp in Fig. 8. 




Fig. 8. The quasiparticle energy bands (scaled by W) relative to the renormalized chemical potential in the t-matrix 
approximation for < fc < 2kF and U/W = 8.0. The states with e < are occupied. 



The resultant density of states, /Oc(e) and Pp(e) in the half-filled case is shown in Fig. 9 for 
U/W = 8.0. 

In the t-matrix approximation, there exists the contribution to the spectral weight also from the 
continuum part, i.e., Pc{k, e) 7^ 0, although its contribution to the density of states, /9c(e)) is small. 
It is seen that Pp{e) has appreciable spectral weight near the Fermi level e ~ 0. We note that for 
U/W » 1, p{e) / over the range of |e| < U/A. 
3.3.3 The Modified T-matrix Approximation 

Finally, we include Umklapp processes by use of Tii{q, z). The modified t-matrix approximation 
leads to S(/c,ie„) for A; > as given by 

T,{k,ien) = -T^^Go{-k + q, -ze„ + iLOi)Tu{q,iLOi). (3.25) 

iuji q 
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Fig. 9. The density of states in the t-matrix approximation for the case of half fiUing and U jW = 8.0. 



In this section, wc will consider only the case of half filling for simplicity. When we calculate 
S(A;,ien) for A; > 0, we take only the account of contributions of the Tii{q,iujl) with q ~ 2kp. 
After straightforward calculations, wc can evaluate e) and then p{k,e). (The first term of 

the r.h.s. of (3.11) is included in the chemical potential shift.) There arc three quasiparticle 
energy bands corresponding to the poles of G^{k,e). These are shown in Fig. 10 for k > and 
U /W = V/W = 8.0. The Umklapp scattering produces the splitting between the upper and lower 
band, between which there exists another energy band. 



UAV=V/W=8.0 




0.5 1 1.5 



Fig. 10. The quasiparticle energy bands (scaled by W) relative to the renormalized chemical potential in the mod- 
ified t-matrix approximation for the case of half filling, A; > and U /W = V/W = 8.0. The states with e < are 
occupied. 



The density of states is shown in Fig. 11 for U/W = 8.0 and a few values of V/W. It is seen that 
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for V ^ 0, p{e) has three-peak-structure. As the Umklapp scattering becomes strong, the upper 
and lower bands approach the Fermi level and the weight in these bands increases, while the weight 
in the middle peak near the Fermi level decreases. 




Fig. 11. The density of states in the half-fiUed case in the modified t-matrix approximation for U /W = 8.0 and 
V/W = 4.0, 8.0. 



The density of states in the case of U = V is shown in Fig. 12. As C/ or F grows up in this 
case, the upper and lower band move away from the Fermi level and three-peak-structure becomes 
prominent. In this case, for U/W » 1, p{e) 7^ over the range of |e| ^ VU /2. 




Fig. 12. The density of states in the modified t-matrix approximation for the case of half filling and U /W = V/W = 
4.0, 20.0. 



As seen, Umklapp processes introduce the global splitting of the density of states, although there 
remains the appriciable spectral weight in the middle of split peaks and we cannot describe the 
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Mott insulator. 

We note that we calculated the self energy by the perturbation theory with the bare Green 
function Gq. In this calculation, the system behaves like the Fermi liquid in that the renormalization 
factor Z ^ 0. We consider that this is due to the lack of the self-consistency. 

But we emphasize that we are interested in the global feature of the density of states, p(e), as a 
function of e in the half-filled case in the presence of the forward and Umklapp scattering. 

§4. Mean Field Approximation 

In the last section, we found that if the Umklapp scattering becomes dominant, there appears a 
pole in the upper half plane in the modified t-matrix. In this section, we examine this instability 
in the mean field approximation. 

We use the cosine band eq. ( 2.2cl|) as the dispersion. Hereafter, for simplicity, we set [/ = and 



consider V only, and examine the ground state properties at absolute zero, T = 0. 
4-1 The G/2— pairing 

The fact that there exists an instability in the particle-particle channel with the total momentum 
\q\ ~ Ikp naturally introduces the following order parameters, 

= 77 2^ < c_fc+G/2iCfcT >, (4.1a) 
k 

V v-^/ 

^L = ^2^ < C-kiCk-cm >, (4.1b) 
k 

where 

k 0<k<G/2 

Aji{Al) is the particle-particle pairing with the total momentum G/2(— G/2). In the half-filled 
case, G/2 is equal to 2kF- 

We extract q = G/2 term from Hjj'', which is decoupled to obtain the mean field Hamiltonian, 

jjMF^ ^ H^ + H,, (4.2a) 

a=R,L 

Ha = Y!'('iaMkai^ka, (4.2b) 
k 

4Re(A^A2) + ^'(e-fc+G/2 + C-k), (4.2c) 

Mkn=\ \,MkL=\ ^TJ' l- (4-2e) 





13 



We note that there exists a self-consistent solution only for the case of 
The quasiparticle dispersion in this case is given by 

^tik) = Caek ± Eg, 



-Ar =Ai. 



Er, 



A^2 + |Ai|2, 



where Cq = 1(— 1) for a = R{L). This is shown in Fig. 13(a). 



(4.3a) 
(4.3b) 



(b) 




G/2 



(C) 



G/2: 



Fig. 13. The quasiparticle energy bands relative to the chemical potential for (a)G/2— pairing state, (b)G'/2— SDW 
state and (c) the coexisting state with both G/2— pairing and G/2— SDW. Full line shows the occupied state. 



It is seen that Ai ^ only for V > 8t, in which case we have 

/. = ^(n-l), (4.4) 

|Ai| = -^^n(2-n), (4.5) 

where n (0 < n < 2) is the filling (the expectation value of the number of electrons) per lattice site. 
(There is another solution for 8t < V < 127rt, although it corresponds to a local maximum of the 
energy.) It turns out from eq. (4.4) and (4.5) that Eg = V/A for V > 8t and the shape of the energy 
band does not change for arbitrary filling n and the Fermi level lies always in the energy gap. Such 
a state is possible, because the weight of the lower band shifts into the upper band as n decreases. 
If we use the linear dispersion, eq. (3.1), nonzero Ai exists for V > 2'kv = 2W . Therefore it is 
seen that this corresponds to the condition that the modified t-matrix has a pole in the upper half 
plane. The essential condition for the stability of this ordered state is + e_k_^Qi2 = 0. 
The total energy per lattice site, Ei , for fixed n is 

El = -^n(n - 2). (4.6) 
o 

We note that the existence of Ai does not lead to the r/-pair condensationEi'0', since 

1 

iV 



Ori-pair = ^ ^{-lYcii^Cil, (4.7) 
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< O^-pair > = -y{AR + At), 
= 0. 



However the following type of coherence does exist, 



^ yn(2-n) (-1)^ 

TT j - jl 

for j/ = j±l,j±3, ■ 



The n-dependences of |Ai|, /x and Ei are shown for V/t = 12 in Fig. 14. 



(4.8) 




Fig. 14. The dependences of |Ai|, chemical potential, /ti, and total energy, E\, (scaled by t) on the electron density 
in G/2— pairing state for V/t = 12 . 



4.2 The G/2-SDW 

Since two electrons with k ^ kpior —kp) are scattered into the final state with k ^ —kpior kp), 
the following particle-hole pairing order parameters are expected to be stabilized: 

^-=^E'<4-G/2.Cfca>, (4.9) 

k 

where a =t, i and J2' f. = J2o<k<G/2- The mean field Hamiltonian by using these order parameters 
is given as follows, 

H^P = Y^H, + H,, (4.10a) 

u 

Ha = E' ^la^kai^ka. H, = -^Re{A^A^), (4.10b) 

k 

^ka ={ I ,Mfe, = I f I , (4.10c) 

\ Ck-G/2a J \ ^ ik-G/2 ) 



15 



where a stands for the spin opposite to a. 

Here we take A| = — A| = A2 so that there exists a self-consistent solution. We write A2 = 
|A2|e*^ where O<0<7r/2. As we shall see later, this phase factor is associated with the amplitude 
of the spin-density-wave (SDW) order parameter. 

The quasiparticle dispersion given by eq. (4.10) is as follows, 

E±{k) = ±E{k) - ^l, (4.11a) 



E{k)^Jel + \/\2?. (4.11b) 



These are shown in Fig. 13(b). In this case the energy gap is 2|A2|. 
We can determine A2 numerically from the self-consistent equation, 



V /-fa-ii-"!) dk , , 



and the chemical potential fi is obtained by using this A2, 



H = sgn(n- l)y'|A2p + At"^ cos^ {^n) . (4.13) 

This solution exists for nc2 < n < 2 — nc2, where 

2 1 r , /47rt., . . 

nc2 = -sin Mtanh( ) . (4.14) 

vr V 

We see from eq. ( [4.12| ) that in the half-filled case (n = 1) this solution exists once V > 0. In the 
limit {V/t = 00), the energy gap is equal to V/2. The nesting property + ^k-G/2 = is again 
crucial for the presence of this order parameter. 
The total energy per lattice site, E2, for fixed n is 

2 /-Kl-ll-nl) 1^ 12 

E2 = / E{k)dk + r-^. (4.15) 

TT Jo y 

In accordance with nonzero A2, the long range order of SDW is induced, but not that of CDW, 
because 

1 

iV 

7(7 

OsDW = ^ ^(-1)V4q<„ (4.17) 

< OcDW > = |rRe(A| + A^), 

= 0, (4.18) 
1 



OcDW^^J2(-^y4aCi., (4.16) 



< OsDW > = |7Re(A| - A|), 



2^cos0. (4.19) 
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We note that the system is degenerate with respect to 9 and in the case 9 = tt/2 the SDW order 
vanishes. 

Recently Wilczek et.al. have introduced the same order parameterEl^ : 

< cl+G/2a^ka >= pure imaginary, 

which corresponds to ^ = 7r/2 in the present notation, i.e., without SDW order. They showed that 
this order (which was named as 'd-Density order') exists in the 2-d extended Hubbard model with 
both on-site and nearest-neighbor Coulomb repulsion, and have discussed the possible relationship 
to the Mott insulators. 

The n-dependences of IA2I, /x and E2 are shown for V/t = 12 in Fig. 15. The important thing is 




Fig. 15. The dependences of IA2I, chemical potential fi and total energy E2 (scaled by t) on the electron density in 
G/2-SDW state for V/t = 12. For n > Uci IA2I / 0. In the half-filled case, = 0. 



that E2 is decreasing function of n , i.e., the compressibility k is negative, k = dn/dfx < 0, except 
for the half-filled case n = 1 at which /i changes discontinuously. Therefore there occurs phase 
separation between a electron-rich and a electron-free phase. 

4.3 The G/2-pairing & G/2-SDW 

It is seen, in fact, that both Ai and A2 can coexist. To see this, we construct the mean-field 
Hamiltonian as follows: 

= Y! ^Mki^k + H„ (4.20a) 

k 

^k = (4t C-fc+G/2i 4-G/2T c-fci)' (^-^Ob) 
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Mk 



Sk-iJ, -Ai 




H. 





-AS 



-A2 






-A2 
Ai 



e--(|Aip + |A2p)-/.. 



We note that, if both Ai and A2 are nonzero, the Cooper-pair also exists in general, 

1 



<Oc> 



N 



E 



(4.20c) 

(4.20d) 

(4.21) 
(4.22) 



-G/2<k<G/2 



^0. 



Hence it is necessary, in principle, to consider the Cooper-pairing term in the mean field Hamilto- 
nian. In the present case of U = 0, however, such a term is given by 



^ c_fe^Cfe| + (h.c.) \ , 



-G/2<k<G/2 

V 



-Gl2<k<G/2 



+ V <Oc> C^G/2^^G/2i + (^-cO^ 

i.e., the term of the total momentum = is associated only with that of = G in Umklapp 
processes. Since the term q = G consists of only one state of two electrons, each of which has wave 
vector G/2, this gives the correction of 0{jj) to the total energy and we can neglect it. 
In the presence of Ai and A2, the quasiparticle energies are given by ±E±{k), 



E±{k) = ^Jel + ^i^ + \Al\^ + \A2\^±2A{k), 
A{k) ^ .Jel\A,\^ + n^{el + \A2n 



(4.23a) 
(4.23b) 



This is shown in Fig. 13(c). 
The self-consistent equations in this case are 



n 



vr 



r 1 1 



A{k) 'E+{k) E^{ky 'E+{k) E_{k) 



Ai = VAi 



■'/'^ dk{ 



1 



r 1 1 



At: \ A{kyE+{k) E_{ky ^E+{k) E^{k)^ j ' 



A, 



VAo — 



dk j 



47r \A{k)^E+{k) E_{k)^ ^ ^E+{k) ^ E_{k) 



1 



+ 



1 



(4.24a) 
(4.24b) 
(4.24c) 
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These are solved numerically. If there is a solution with Ai,A2 7^ 0, the correlation of the 
Cooper-pair, eq. (4.22), is obtained as follows: 

r'^/^ dk 1 , 1 1 , 



< 6c >= -^|Ai||A2|sini 



r7\ 



2TrA{k)'E+{K) E^{ky ^^'^^^ 
We note that the phase factor, 6, associated with the amplitude of SDW order parameter, eq. 
(4.19), appears also in the correlation of the Cooper-pair. 
The solution with Ai, A2 7^ exists for rid < n <2 — rid ifV> 8t, where 



nd = l- 



\ 



1 



(4.26) 



In the limit of F ^ 8t, rid and Ai 0. It is seen that for a half-filled band n = 1 {jj, = 0), 
this solution exists only in the limit t = 0. 
The total energy per lattice site, E^, for fixed n is 

/o 

The n-dependences of |Ai|, IA2I, |Ac| = V\< Oc > |, and E3 are shown in Fig. 16 for V/t = 12 
and 6 = 7r/2. In the case of Ai, A2 7^ 0, the compressibility k becomes positive in contrast to the 
case of Ai = and A2 7^ 0. 



/■7r/2 fiu 2 
E3 = - - [E+{k) + E_{k)] + -(|Ai|2 + IA2P) + ^i{n - 1). 
Jn TT V 



(4.27) 



-1 



®- s- -s n-a-" ' 

t-Al ^-:^^^^^^^^4A-l- 

■/ ^ 

cli . i . 

0.5 0.6 0.7 0.8 0.9 

n 



Fig. 16. The dependences of |Ai|, IA2I, |Ac|, chemical potential n and total energy E3 (scaled by t) on the electron 
density in the coexisting state with both G/2— pairing and SDW for V/t = 12 and 6 = -k/2 . These two states 
coexist for rici < n < 1. 



4.4 Phase Diagram 

The state with the lowest energy has been determined in the V/t — n plane for the region of 
< n < 1. The result is shown in Fig. 17. 
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Fig. 17. Phase diagram in the plane of V jt and n for C/ = 0. PS stands for phase separation. Line AB stands for 
Wc. Line BCD stands for the boundary between the normal state and the coexisting state with both G/2— pairing 
and G/2— SDW. For detailed explanations about tIc, Wci and nc2, see the text. 



To begin with, we consider the case V jt > 8. We compare G/2— pairing state with the normal 
state. By noting that the total energy of the normal state Eq is equal to — ^ sin(^n) and comparing 
this with El, eq. ( [4. 61) , we see that the boundary between the normal state and the G/2— pairing 
state is of the first order as V is varied and is given as follows, 

< n < ric for 8t < y < Vc, 

(4.28) 

< n < 1 for y > Vc, 

where Vc/t is equal to 32/7r(~ 10.2) and ric is determined as the solution of the following equation, 

V At IT 

— ricK - 2) = sin(-nc). (4.29) 

o vr z 

Next we discuss G/2— SDW state. The transition from the normal state to this state is of the 
second order as a function of V. As discussed in the last section, this state is unstable against the 
phase separation (except for the case of half filling). However, it is energetically favorable compared 
to the normal state, E2 < Eq for < n < 1, where E2 and is given by eq. ( [4.15 ) and eq. 
( [4.14 ), respectively. On the other hand, this state competes with G/2— pairing state for V > Vd 
and nc2 < n < rid, where Vd/t ~ 13.0 and rid is given by eq. ( [4.26| ). (Vd is determined as the 
solution of Ud = nc2-) It is seen that in this region G/2— pairing state is stable, Ei < E2- 

In the coexisting state with both G/2— pairing and G/2— SDW for V/t > 8 and Ud < n < 1, the 
total energy, E^, eq. ( 4.27[ ), is always lower than that of the state with only one order parameter 
among them, i.e., E3 < Ei,E2- On the other hand, for V/t ~ 8, the energy of this state is higher 
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than that of the normal state, Eq < E^, i.e., the transition from the normal state to this coexisting 
state is of the first order as a function of V. The boundary between the normal state and this 
coexisting state is determined numerically as the solution of Eq = E^. It is seen that E^ < Eq, 
independent of n, for V/t ^ 9.23. Especially in this region, the ground state is given as follows: 

n = 1 G/2- SDW, 

na < n < 1 G/2 - pairing & G/2 - SDW, (4.30) 
< n < rici G/2 — pairing, 

On the other hand, for < F/i < 8, the ground state is given as follows: 

n = 1 G/2- SDW, 

' (4.31) 
< n < nc2 normal, 

and for nc2 < n < 1 there occurs phase separation. 

§5. Conclusion and Discussion 

Wc have studied the effects of the Umklapp scattering in the ID Hubbard model by treating 
the strength of normal (U) and Umklapp (V) interaction as independent parameters, keeping the 
bandwidth (W) finite. 

We calculated the density of states by use of the perturbation theory with the linear dispersion 
with a special emphasis on both the forward and Umklapp scattering. 

For V = 0, we found that the t-matrix (the summation of the particle-particle ladder diagram 
with respect to U) has one pole on the real axis. The density of states, p{e), in the t-matrix 
approximation, has appreciable spectral weight and no energy gap at the Fermi level in the half- 
filled case. 

In the presence of V, we found that the modified t-matrix (including V term) has two poles on 
the real axis for < V < U . In the resultant p{e) in the modified t-matrix approximation, the 
global splitting is introduced and three-peak-structiirc is observed in the half-filled case, although 
the appreciable spectral weight remains in the middle of split peaks near the Fermi level. As the 
Umklapp scattering becomes strong, this structure becomes prominent. However, we have not 
described the Mott insulator. 

On the other hand, ii V > U + 2W, a pole of the modified t-matrix appears in the upper half 
plane. Therefore in §4, first of all wc studied this instability in the mean field approximation for 
U = and found that in the strong coupling regime, G/2— pairing state becomes stable. Next 
we examined another instability associated with Umklapp processes and found that there appears 
G/2— SDW state. This state, however, is mctastable due to the negative compressibility except for 
a half-filled band. Finally, we found that the above two order parameters can coexist. 

We constructed a phase diagram in the full V/t — n plane for these possible states. In the strong 
coupling regime V/t > 8, the ground state is the G/2— SDW state at just half filling n = 1, the 
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coexisting state with both G/2— pairing and G/2— SDW near half fiUing and the G/2— pairing state 
at low filling. In this coexisting state near half filling, the correlation of the Cooper-pair appears 
and is connected with the SDW order parameter through one phase angle , 6, i.e., < Oc >oc sin 6*, 
eq. ( 4.25| ), and < OsDW >oc cos9, eq. (4.19), through the existence of the G/2— pairing. Such 



close relationship between the SDW and Cooper-pair has been indicated in the 2D t — J model 
with the S0(5) symmetrylll) , which unifies antiferromagnetism with d-wave superconductivity. The 
application of the present theory to such 2D case is the future problem. 
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